Comparison of Three Motion Capture-Based Algorithms for Spatiotemporal Gait Characteristics: How Do Algorithms Affect Accuracy and Precision of Clinical Outcomes?

Gait assessment is of interest to clinicians and researchers because it provides information about patients’ functional mobility. Optoelectronic camera-based systems with gait event detection algorithms are considered the gold standard for gait assessment. Yet, the choice of the algorithm used to process data and extract the desired parameters from those detected gait events has an impact on the validity and reliability of the gait parameters computed. There are multiple techniques documented in the literature for computing gait events, including the analysis of the minimal position of the heel and toe markers, the computation of the relative distance between sacrum and foot markers, and the assessment of the smallest distance between the heel and toe markers. Validation studies conducted on these algorithms report variations in accuracy. Yet, these studies were conducted in different conditions, at varying gait velocities, and on different populations. The purpose of this study is to compare accuracy, precision, and robustness of three algorithms using motion capture data obtained from 25 healthy persons and 21 psoriatic arthritic patients walking at three distinct speeds on an instrumented treadmill. Errors in gait events recognition (heel strike—HS and toe-off—TO) and their impact on gait metrics (stance phase and stride length) are reported and compared to ground reaction force events measured with force plates. Over the 9114 collected steps across all walking speeds, more than 99% of gait events were recognized by all algorithms. On average, HS events were detected within 1.2 ms of the reference for two algorithms, while the third one detected HS late, with an average detection error of 40.7 ms. Yet, significant variations in accuracy were noted with gait speed; the performance decreased for all algorithms at slow speed. TO events were identified early by all algorithms, with an average error ranging from 16.0 to 100.0 ms. These gait events errors lead to 2–15% inaccuracies in stance phase assessment, while the impact on stride length remains below 0.3 cm. Overall, the algorithm based on the relative distance between the sacral and foot markers stood out for its accuracy, precision, and robustness at all walking speeds.


Introduction
Functional mobility is the ability to move independently and safely in all types of contexts to execute everyday activities, such as walking, bending, rising from a chair, performing housework, and dressing up. Loss of functional mobility, whether owing to aging or disease, has a detrimental effect on the performance of these activities and can impact the quality of life of the individual [1]. Clinicians and researchers are, therefore, Thus, the scientific literature reveals that algorithms used for gait event detection perform differently, which in turn can impact the accuracy of the gait performance parameters extracted. Yet, several studies have used optoelectronic-based algorithms interchangeably either to validate another detection method or to answer a clinical question [15,[18][19][20][21][22]. This is concerning as gait event detection algorithm robustness and accuracy appear to vary with population and gait velocity. Moreover, current investigations on clinical populations were conducted under varying conditions and on a small number of individuals, making it difficult to combine the results and obtain a clear picture. The therapeutic implications of these limitations (e.g., stance phase accuracy and step length) are not well understood either. The latter is essential to accurately assess abnormal gait (usually slower gait), as well as the progression of a disorder or the effect of an intervention.
To the best of our knowledge, no work has examined and compared gait event algorithms based on optoelectronic markers at different walking speeds to evaluate the systematic inaccuracy in gait event detection timings and their impacts on gait metrics used to characterize gait. Specifically, how much variation in gait event detection is attributed to these regularly used algorithms? What effect does this variation have on clinical metrics?
The objectives of the current study were, thus, to evaluate the accuracy of three MoCapbased algorithms at three distinct walking speeds to measure gait events (HS and TO) and assess their impact on gait metrics (stance duration, stance phase, and stride length), computed in patients with psoriatic arthritis (PsA) and healthy adults. This research aims to contribute to a greater understanding of actual variances, as well as a determination of which algorithm is the most robust for adoption as the official gold standard.

Materials and Methods
This study was approved by the CIUSSS de l'Estrie-CHUS ethics committee and conducted at the Centre de Recherche sur le Vieillissement (CdRV), within the Université de Sherbrooke, Canada. All participants provided written informed consent.

Participants
A total of 46 participants were enrolled in the study, including 25 healthy adults (HP) (15 females, mean age of 31.1 ± 10.2 years, mean body mass index (BMI) of 24.6 ± 4.3) and 21 persons living with foot impairments secondary to psoriatic arthritis (PsA) (16 females, mean age of 53.9 ± 8.9 years, BMI of 29.3 ± 4.5, mean disease duration of 11.5 ± 10.2 years). For HP, participants had to be over 18 years old and in good general health to qualify. People with lower-limb pain or any musculoskeletal, rheumatological, or neurological disease that could affect normal gait patterns were excluded from the study. Participants with PsA were recruited sequentially from the rheumatology outpatient clinics at the Hotel Dieu University Hospital CHU of Sherbrooke (CHUS). The inclusion criteria were as follows: aged between 20 and 70 years, a confirmed diagnosis of PsA by a trained rheumatologist, moderate to severe and recurrent foot pain, and stable medication over 3 months preceding the recruitment. Patients with diabetes, neurological disease, or any musculoskeletal ailment that could impact normal gait patterns were excluded.

Measurement Systems
Participants' gait was evaluated as they walked on a split-belt instrumented treadmill (AMTI, Watertown, MA, USA) equipped with two force plates (sampling rate = 1000 Hz) for 2 min, at slow, normal, and fast speed (0.45, 1.12, and 1.6 m/s) [23]. The order of the test condition was balanced with a Latin square design. Participants were instrumented with 16 passive markers following the requirements for the use of OptiTrack's standard lowerbody model [24]. Specifically, the chosen model includes markers on the lateral prominence of the ankles malleolus bones, on the dorsal face of the second metatarsi, and on the calcaneus bone of the feet. The 3D position of the markers was captured at 100 Hz using an Optitrack Motion Capture (MoCap) system (NaturalPoint, Corvallis, OR, USA) composed of eight Prime 13 W cameras, carefully positioned and calibrated according to the manu- facturer's instructions to achieve the claimed 0.30 mm tracking accuracy over the desired volume of acquisition ( Figure 1) [25].
test condition was balanced with a Latin square design. Participants were instrumented with 16 passive markers following the requirements for the use of OptiTrack's standard lower-body model [24]. Specifically, the chosen model includes markers on the lateral prominence of the ankles malleolus bones, on the dorsal face of the second metatarsi, and on the calcaneus bone of the feet. The 3D position of the markers was captured at 100 Hz using an Optitrack Motion Capture (MoCap) system (NaturalPoint, Corvallis, OR, USA) composed of eight Prime 13 W cameras, carefully positioned and calibrated according to the manufacturer's instructions to achieve the claimed 0.30 mm tracking accuracy over the desired volume of acquisition ( Figure 1) [25]. Figure 1. Study overview. Participants performed three 2 min treadmill walking trials at slow, normal, and fast walking speeds on a split-belt instrumented treadmill. Participants were instrumented with 16 passive markers on the lower body. Gait was simultaneously captured by both the treadmill force plates and the motion capture system. Gait events were then assessed using three different motion capture algorithms.
The treadmill and the MoCap system were synchronized using an eSync hub, ensuring a synchronization accuracy within 5 ms. Calibration and data acquisition were performed with Motive: Body (Version 2.1.1). Missing data within optoelectronic markers' position signals were interpolated using a shape-preserving piecewise cubic interpolation approach.

Algorithms
This section describes in detail the algorithms used to detect gait events (HS and TO) using force plates, our gold standard, and kinematic data. Recall that HS is defined as the initial contact of the foot to the ground, and TO represents the instant corresponding to the last contact of the foot with the ground. Hence, heel strike is referred to as the beginning of the gait cycle, and toe-off marks the end of the stance phase [7].

Gold Standard: Force Plates (GS)
Gait event detection based on force plate data is considered the gold standard [14]. HS and TO represent initiation and termination instants of the foot contact to the ground, which corresponds to the change in detected force levels. Identification of HS and TO (i.e., the beginning of the rising edge and the end of the falling edge of the force signal, respectively) is often performed using a fixed (5 N [26], 10 N [27][28][29], 20 N [14,22], or adaptative percentage of the mass [30]) threshold. In the current study, these instants were detected Figure 1. Study overview. Participants performed three 2 min treadmill walking trials at slow, normal, and fast walking speeds on a split-belt instrumented treadmill. Participants were instrumented with 16 passive markers on the lower body. Gait was simultaneously captured by both the treadmill force plates and the motion capture system. Gait events were then assessed using three different motion capture algorithms.
The treadmill and the MoCap system were synchronized using an eSync hub, ensuring a synchronization accuracy within 5 ms. Calibration and data acquisition were performed with Motive: Body (Version 2.1.1). Missing data within optoelectronic markers' position signals were interpolated using a shape-preserving piecewise cubic interpolation approach.

Algorithms
This section describes in detail the algorithms used to detect gait events (HS and TO) using force plates, our gold standard, and kinematic data. Recall that HS is defined as the initial contact of the foot to the ground, and TO represents the instant corresponding to the last contact of the foot with the ground. Hence, heel strike is referred to as the beginning of the gait cycle, and toe-off marks the end of the stance phase [7].

Gold Standard: Force Plates (GS)
Gait event detection based on force plate data is considered the gold standard [14]. HS and TO represent initiation and termination instants of the foot contact to the ground, which corresponds to the change in detected force levels. Identification of HS and TO (i.e., the beginning of the rising edge and the end of the falling edge of the force signal, respectively) is often performed using a fixed (5 N [26], 10 N [27][28][29], 20 N [14,22], or adaptative percentage of the mass [30]) threshold. In the current study, these instants were detected using a plateau edge enhancement technique. The force plate signal was first filtered using a third-order low-pass Butterworth filter with a cutoff frequency of 50 Hz. The filtered signal was then translated in order to position the low-level plateau (i.e., foot not in contact with the ground) away from zero. The modified signal was then multiplied by an exponential function (decreasing/increasing exponential to enhance rising/falling edge, respectively). HS and TO were then identified using a peak detection approach on the enhanced signal (see Figure 2i). using a plateau edge enhancement technique. The force plate signal was first filtered using a third-order low-pass Butterworth filter with a cutoff frequency of 50 Hz. The filtered signal was then translated in order to position the low-level plateau (i.e., foot not in contact with the ground) away from zero. The modified signal was then multiplied by an exponential function (decreasing/increasing exponential to enhance rising/falling edge, respectively). HS and TO were then identified using a peak detection approach on the enhanced signal (see Figure 2i).

Minimal Position of the Heel and Toe Markers (MPHT)
The first MoCap-based algorithm was a modified version of Hreljac and Marshall's algorithm [26]. It is based on an analysis of the minimal vertical position of the heel (calcaneus) and toe (second metatarsal) markers. A third-order low-pass Butterworth filter with a cutoff frequency of 15 Hz was first applied to the signals. HS was determined as the heel minimum vertical position. Similarly, terminal contact is identified as the toe minimum vertical position, corresponding to the propulsion and, thus, to the initiation of vertical displacement (see Figure 2ii).

Relative Distance between Sacral Marker and Foot Markers (RDSF)
The second MoCap-based algorithm considers the relative distance between a foot marker and the sacral marker using the frontal axis of the heel and toe markers' position [15]. The heel and toe signals were filtered with a third-order low-pass Butterworth filter, using a cutoff frequency of 15 Hz. Figure 2iii shows the sinusoidal curve resulting from the analysis of the distance between the heel and toe markers with the sacral marker projected onto the estimated walking direction. Peaks and valleys with amplitude over 30% of the mean signal amplitude were determined as HS and TO events, respectively. and toe (red) markers (MPHT); (iii) the resultant frontal coordinate formed by the subtraction of the frontal coordinate of the sacral marker from the frontal coordinate of the heel marker (blue) and resultant frontal coordinate formed by the subtraction of the frontal coordinate of the sacral marker from the frontal coordinate of the toe marker (red) (RDSF); (iv) the resulting signal formed by the subtraction of the vertical coordinate of the heel and toe markers (MDHT). In all subfigures, blue squares represent detected heel strikes, and red circles represent identified toe-offs.

Minimal Position of the Heel and Toe Markers (MPHT)
The first MoCap-based algorithm was a modified version of Hreljac and Marshall's algorithm [26]. It is based on an analysis of the minimal vertical position of the heel (calcaneus) and toe (second metatarsal) markers. A third-order low-pass Butterworth filter with a cutoff frequency of 15 Hz was first applied to the signals. HS was determined as the heel minimum vertical position. Similarly, terminal contact is identified as the toe minimum vertical position, corresponding to the propulsion and, thus, to the initiation of vertical displacement (see Figure 2ii).

Relative Distance between Sacral Marker and Foot Markers (RDSF)
The second MoCap-based algorithm considers the relative distance between a foot marker and the sacral marker using the frontal axis of the heel and toe markers' position [15]. The heel and toe signals were filtered with a third-order low-pass Butterworth filter, using a cutoff frequency of 15 Hz. Figure 2iii shows the sinusoidal curve resulting from the analysis of the distance between the heel and toe markers with the sacral marker projected onto the estimated walking direction. Peaks and valleys with amplitude over 30% of the mean signal amplitude were determined as HS and TO events, respectively.

Minimal Distance between the Heel and Toe Markers (MDHT)
The last MoCap-based algorithm finds its grounds in the theoretical definition of gait cycle and events [31,32]. Heel and toe positions were first filtered using a thirdorder low-pass Butterworth filter with a cutoff frequency of 15 Hz. The relative vertical distance between heel and toe markers was then analyzed to detect HS (heel-toe). Indeed, HS corresponds to the minimum vertical relative distance and is detected with Matlab's findpeaks function. TO was associated with the final burst of propulsion of the gait cycle. It was, thus, detected at the peak rate of change in relative distance between two consecutive Sensors 2023, 23, 2209 6 of 13 heel strikes. Figure 2iv illustrates an example of a resulting signal and its associated HS and TO.

Data Analysis
For each trial, the first and last 20 s were discarded to focus on steady-state walking. The remaining 80 s were then analyzed to identify heel strikes (HS) and toe-offs (TO) from the force plate signals (gold standard) and the motion-capture system (three different algorithms). The four methods for gait event detection were implemented and processed automatically in Matlab R2020b.
Accuracy was assessed using the mean and standard deviation of HS and TO error in terms of detection time for each walking speed and algorithm. Robustness of the algorithms was evaluated on the basis of their capacity to detect gait events and extract spatiotemporal gait parameters. Stance duration refers to the time elapsed between consecutive HS and TO on a specific side. Stance phase is a normalized version of stance phase duration over the complete gait cycle duration, expressed as a percentage. Stance duration is, thus, useful to assess measurement accuracy, while stance phase, as a percentage, refers to its clinical assessment. Stride length corresponds to the distance the foot traveled between two consecutive heel strikes. It was appraised as the heel marker's displacement between two consecutive HS events to which the product of the treadmill speed by the time elapsed between the two HS events was added. The mean value, errors, and bias of those parameters were used to evaluate the algorithms' accuracy, while the precision was evaluated using the standard deviation. The effects of the algorithm and the walking speed on gait event detection and clinical parameter accuracy were then assessed using a two-way ANOVA. Factor interactions were explored using Scheffé's test. Statistical significance was considered with an α = 0.05. All statistical analyses were performed using SPSS Statistics 26.
The reliability of algorithms was assessed through an evaluation of the limits of agreement between the specific algorithm and the gold standard using a Bland and Altman approach. Clinical impact was evaluated in terms of reliability of clinical parameters (stance duration and stride length), assessed using the standard error measurement (SEM) and minimal detectable change (MDC) approaches [33]: where r is the reliability coefficient SEM and MDC were calculated considering each algorithm compared to the gold standard to assess the reliability of the algorithms. Recall that the reliability coefficient varies between 0 and 1, where a high value is associated with a high reproducibility. There is no real consensus regarding "cutoff" values for the interpretation of reliability coefficients. However, it has been suggested that values <0.59 reflect "poor reproducibility", values of 0.60-0.79 reflect "moderate reproducibility", and values >0.80 reflect "high reproducibility" [34].

Gait Events Detection
A total of 138 trials were collected during this study. However, one healthy participant's slow trial was excluded from the analysis because of irregularities in the force plates signals. For PsA, participants were not able to reach the desired velocity or to maintain it for 2 min on 34 occasions. Thus, a total of 103 trials (74 for HP: 24 slow speed, 25 normal speed, 25 fast speed; 29 for PsA: 15 slow speed, 11 normal speed, three fast speed) including 9114 steps were analyzed using all four algorithms (GS, MPHT, RDSF, and MDHT). All three MoCap-based algorithms detected more than 99% of gait events (HS and TO).
Distribution of the timing difference in HS and TO detected with the three MoCap-based algorithms compared to GS is shown in Figure 3 at slow, normal, and fast walking speeds. A positive value indicates a delay in event detection for the optoelectronic method in comparison to the GS.
A total of 138 trials were collected during this study. However, one healthy participant's slow trial was excluded from the analysis because of irregularities in the force plates signals. For PsA, participants were not able to reach the desired velocity or to maintain it for 2 min on 34 occasions. Thus, a total of 103 trials (74 for HP: 24 slow speed, 25 normal speed, 25 fast speed; 29 for PsA: 15 slow speed, 11 normal speed, three fast speed) including 9114 steps were analyzed using all four algorithms (GS, MPHT, RDSF, and MDHT). All three MoCap-based algorithms detected more than 99% of gait events (HS and TO). Distribution of the timing difference in HS and TO detected with the three MoCap-based algorithms compared to GS is shown in Figure 3 at slow, normal, and fast walking speeds. A positive value indicates a delay in event detection for the optoelectronic method in comparison to the GS. The mean and standard deviations (SD) in detection times errors for all algorithms and walking speeds are reported in Table 1. Again, a positive value indicates that the optoelectronic method detected the event later than GS. On average, MPHT detected HS later and TO detected HS prematurely for all walking speeds and both populations. Average HS was detected on time with both RDSF and MDHT algorithms. Yet, RDSF detected TO closer to the reference value than MDHT, which detected it prematurely. It is also interesting to look at the distribution of HS and TO accuracy in terms of frame differences between the detected event and the gold standard. Overall, 59.6% of HS events were detected within one frame using MDHT algorithms, and more than 90% were identified within three frames. While similar results were achieved with RDSF (57.1% within one frame, 93.2% within three frames), MPHT detected only 1.5% of HS events within one frame, and the 90% detection level was reached within five frames. RDSF performance The mean and standard deviations (SD) in detection times errors for all algorithms and walking speeds are reported in Table 1. Again, a positive value indicates that the optoelectronic method detected the event later than GS. On average, MPHT detected HS later and TO detected HS prematurely for all walking speeds and both populations. Average HS was detected on time with both RDSF and MDHT algorithms. Yet, RDSF detected TO closer to the reference value than MDHT, which detected it prematurely. It is also interesting to look at the distribution of HS and TO accuracy in terms of frame differences between the detected event and the gold standard. Overall, 59.6% of HS events were detected within one frame using MDHT algorithms, and more than 90% were identified within three frames. While similar results were achieved with RDSF (57.1% within one frame, 93.2% within three frames), MPHT detected only 1.5% of HS events within one frame, and the 90% detection level was reached within five frames. RDSF performance with regard to the TO accuracy distribution was very similar to HS (59.7% within one frame, 86.3% within three frames). However, both MPHT and MDHT identified more than 85% of TO events with a difference greater than five frames from the gold standard.
Accuracy variations on TO were then assessed using the same approach. The tests exposed a statistically significant overall difference (p < 0.001), with a significant effect of both factors (p < 0.001). Post hoc tests revealed statistically significant differences for all speed combinations (p < 0.001 for all three combinations; mean differences: slow vs. normal: 34 ms, slow vs. fast: 47 ms, normal vs. fast: 13 ms). Similarly, all three comparisons in algorithms were statistically significant (p < 0.001 for all three; mean differences: MPHT vs. RDSF: 80 ms, MPHT vs. MDHT: 10 ms, MDHT vs. RDSF: 71 ms). Variation in HS accuracy, assessed with a two-way ANOVA on combined HP and PsA data, revealed a statistically significant variation in accuracy (p < 0.001), with a significant effect of both velocity (p < 0.001) and algorithm (p < 0.001). Specifically, post hoc analyses exposed that HS accuracy was significantly different at slow speed (slow vs. normal: p < 0.001; slow vs. fast: p < 0.001). Yet, no significant difference was noted between normal and fast gait as far as HS accuracy is concerned (p = 0.221). However, power analysis on this specific comparison exposed a 39% power. Regarding algorithms, post hoc analyses revealed a significant difference between MPHT and both RDSF and MDHT (p < 0.001), but not between RDSF and MDHT (p = 0.762).

Clinical Impact
Stance duration, stance phase percentage, and stride length were then assessed using GS and all three MoCap-based algorithms to better evaluate the impact of gait events detection accuracy on clinical outcomes. Bland-Altman (BA) analysis on stance duration is presented in Figure 4A. All three MoCap-based approaches revealed a bias in stance duration estimation compared to GS. Yet, the importance of that bias varies with the approach, as does the dispersion, leading to different 95% confidence intervals (MPHT: −282.2 to 19.7 ms; RDSF: −94.3 to 68.7 ms; MDHT: −185.0 to 15.5 ms). Details on stance duration values and errors are listed in Table 2. Two-way ANOVA performed on stance duration revealed a statistically significant variation (p < 0.001). Sheffé's test exposed a significant impact of all gait speed combinations and algorithm comparisons (p < 0.001), although the importance of that impact varied (mean difference in accuracy between gait speeds: slow vs. normal: 400 ms, slow vs. fast: 495 ms, normalvs. fast: 95 ms; mean difference in accuracy between algorithms: RDSF vs. MPHT: 267 ms, RDSF vs. MDHT: 67 ms, MDHT vs. MPHT:200 ms). In terms of stance phase percentage, this led to errors varying between 2% and 15%.
The reliability coefficient (r), SEM, and MDC for both stance duration and stride length clinical parameters based on the impact of the algorithm accuracy are presented in Table 3. Recall that the SEM and MDC based on algorithm accuracy consider the correlation between the MoCap algorithms and GS parameters, thus representing the algorithms' reliability. BA analysis on stride length ( Figure 4B) showed high agreement between the three MoCap-based algorithms and the gold-standard force plate system for all walking speeds. Indeed, bias was lower than 0.1 mm, and confidence intervals were relatively constant (healthy: MPHT: −5.0 to 5.0 cm; RDSF: −4.0 to 4.0 cm; MDHT: −3.9 to 3.9 cm). Two-way ANOVA based on algorithm and walking speed revealed a statistically significant variation (p < 0.001), with a significant impact of walking speed (p < 0.001), but not of the algorithm (p = 0.175). Post hoc analyses exposed a statistically significant difference between all combinations of gait speeds, albeit to a varying extent (p < 0.001 for all pairs; mean difference: slow vs. normal: −458 ms, slow vs. fast: −741 ms, normal vs. fast: −283 ms).

Discussion
In this study, three optoelectronic motion capture-based algorithms were compared to assess their ability to identify gait events during treadmill walking in healthy adults and PsA patients with foot abnormalities. When compared to a force plate gold standard, all MoCap-based algorithms demonstrated great validity, identifying 99% of HS and TO occurrences. However, all three algorithms detected TO prematurely, with a considerably smaller difference observed with the RDSF algorithm, and a general tendency for that error to diminish with increasing walking speed. On average, HS was detected within 2 ms of actual events using RDSF or MDHT, while mean error on HS using MPHT was, on average, 41 ms. Thus, the error observed on clinical parameters was mainly induced by the accuracy of TO events. It is interesting to note that the results obtained for RDSF gait event accuracy were consistent with previously published work by Zeni, Richards, and Higginson which mentions a timing accuracy within 16.7 ms for 94% of gait events [15]. Furthermore, analysis of variance on the accuracy of both HS and TO revealed a significant impact of velocity and algorithm.
Regarding clinical parameters, all three algorithms showed high agreement with the gold standard for stride length, with a global bias close to 0 and overall confidence intervals varying from ±4 cm to ±5 cm As far as stance duration is concerned, agreement varied amongst algorithms; both MPHT and MDHT introduced a bias and had larger confidence intervals than RDSF. With a mean stance phase error within 2% and a mean stride length error below 0.3 cm across gait speed and populations, the RDSF algorithm stood out. Evaluation of the minimal detectable change considers the variability within the data to estimate the minimal true change one can expect to capture with a specific system or algorithm. The MDC for stance phase duration was 40 ms for RDSF and MDHT, while it reached 50 ms for MPHT. In other words, changes in stance duration above 40 ms detected with RDSF or MDHT can be considered as a true change. As such, changes observed in stance duration with either gait speed variation or pathological condition can be assessed properly with these algorithms, according to the current study data. For stride length, although global mean errors were close to zero for all three algorithms, MDC values ranged between 9 cm and 17 cm. This may be explained by the varying importance of the errors to the actual measure. In usual clinical gait assessment cases, gait velocity is closer to the slow or normal conditions of this study. As such, relevant MDC can be considered closer to 4 cm for both RDSF and MDHT. Throughout the present study, the RDSF algorithm, therefore, stood out for its simplicity and reliability in terms of HS and TO detection times, leading to a more accurate estimation of clinical parameters. This algorithm is less sensitive to noise or irregularities in the signal and is simpler to implement. Moreover, accuracy of the RDSF algorithm is better and steadier across all walking speeds, and its MDCs were, overall, the best. Thus, the RDSF algorithm is perceived to be the most efficient algorithm to use with optoelectronic system for gait assessment.
However, as mentioned in the introduction, MoCap-based algorithms have been used as a reference to validate other technologies or approaches in the literature [18,20,22]. As such, a spinoff of this study would be to define an early acceptable level of variation for new technology validation. Indeed, the 100 ms difference in HS detection error across the three MoCap-based algorithm gives a sense of the level of accuracy that researchers are currently dealing with and accepting. These guidelines could, thus, be used for early assessment of a new algorithm's potential.
Some limitations apply to the current study. Even the slowest walking speed was too fast for some participants to complete the 2 min trial, and only three PsA participants were able to complete the 2 min trial at the fastest walking speed. Future studies should, therefore, consider the validity and reliability of RDSF using less standardized environments and include other kinds of gait-impaired participants.

Conclusions
The present study highlighted the differences in accuracy and robustness of three algorithms using an optoelectronic system compared to a gold-standard force plate system, in healthy and PsA participants. The differences highlighted in both gait event detection and clinical parameter accuracy raise a flag regarding the importance of the algorithm when considering the optoelectronic system as a reference. Among the tested algorithms, analysis of the relative distance between sacrum and foot (RDSF) is considered the most accurate and robust. Indeed, HS and TO events extracted with RDSF were more accurate, leading to smaller errors in clinical parameters including stance phase and stride length.  Informed Consent Statement: Informed consent was obtained from all participants involved in the study.

Data Availability Statement:
The data presented in this study are available on request from the author responsible for data curation. The data are not publicly available dur to ethical restrictions.